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ABSTRACT 

New JCMT HARP CO 3-2 observations of the W5 star forming complex are 
presented, totaling an area of 12000 arcmin^ with sensitivity better than 0.1 K 
per 0.4 km s~^ channel. We discovered 55 CO outflow candidates, of which 40 are 
associated with W5 and 15 are more distant than the Perseus arm. Most of the outflows 
are located on the periphery of the W5 HH region. However, two outflow clusters 
are > 5 pc from the ionization fronts, indicating that their driving protostars formed 
without directly being triggered by the O-stars in W5. We compare the derived outflow 
properties to those in Perseus and find that the total W5 outflow mass is surprisingly 
low given the cloud masses. The outflow mass deficiency in the more massive W5 cloud 
(M(H2) ~ 5 X 10^ M0) can be explained if ionizing radiation dissociates molecules as 
they break out of their host cloud cores. Although CO J=3-2 is a good outflow tracer, 
it is likely to be a poor mass tracer because of sub-thermal line excitation and high 
opacity, which may also contribute to the outflow mass discrepancy. It is unlikely 
that outflows could provide the observed turbulent energy in the W5 molecular clouds 
even accounting for undetected outflow material. Many cometary globules have been 
observed with velocity gradients from head to tail, displaying strong interaction with 
the W5 HH region and exhibiting signs of triggered or revealed star formation in their 
heads. Because it is observed face-on, W5 is an excellent region to study feedback 
effects, both positive and negative, of massive stars on star formation. 

Key words: ISM: jets and outflows — ISM: kinematics and dynamics — ISM: indi- 
vidual: W5 — stars: formation 



1 INTRODUCTION 

Galactic-scale shocks such as spiral density waves promote 
the formation of giant molecular clouds (GMCs) where mas- 
sive stars, star clusters, and OB associations form. The mas- 
sive stars in such groups can either disrupt the surrounding 
medium or promote further star formation. While ionizing 
and soft UV radiation, stellar winds, and eventually super- 
nova explosions destroy clouds in the immediate vicinity of 
massive stars, as the resulting bubbles age and decelerate, 
they can also trigger further star formation. In the "collect 
and collapse" scenario (e.g. [Elmegreen Lada 1977), gas 
swept-up by expanding bubbles can collapse into new star- 
forming clouds. In the "radiation-driven implosion" model 



(Bertoldi McKee 1990 Klein et al. 1983), pre-existing 



clouds may be compressed by photo-ablation pressure or by 
the increased pressure as they are overrun by an expand- 
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ing shell. In some circumstances, forming stars are simply 
exposed as low-density gas is removed by winds and radi- 
ation from massive stars. These processes may play signifi- 
cant roles in determining the efficiency of star formation in 
clustered environments ( |Elmegreen||1998 ). 

Feedback from low mass stars may also control the 



shape of the stellar initial mass function in clusters ( Adams 
Fatuzzo]|1996 Peters et al.||2010i . Low mass young stars 
generate high velocity, collimated outflows that contribute 
to the turbulent support of a gas clump, preventing the 
clump from forming stars long enough that it is eventually 
blown away by massive star feedback. It is therefore im- 
portant to understand the strength of low-mass protostellar 
feedback relative to other feedback mechanisms. 

Outflows are a ubiquitous indicator of the presence of 



ongoing star formation (Reipurth & Bally 2001). CO out- 



flows are an indicator of ongoing embedded star formation 
at a younger stage than optical outflows because shielding 
from the interstellar radiation field is required for CO to 
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survive. Although Herbig- Haro shocks and H2 knots reveal 
the locations of the highest- velocity segments of these out- 
flows, CO has typically been thought of as a "calorimeter" 
measuring the majority of the mass and momentum ejected 
from protostars or swept up by the ejecta ( Bachiller|1996 ). 

The W5 star forming complex in the outer galaxy is a 
prime location to study massive star formation and trigger- 
ing. The bright-rimmed clouds in W5 have been recognized 
as good candidates for ongoing triggering by a number of 
groups ([Lefloch et ar]|1997[ [Thompson et al.|[2004l [Karr k 



Martin 2003). The clustering properties were analyzed by 
Koenig et al. ( |2008| ) using Spitzer infrared data, and a num- 



ber of significant clusters were discovered. The whole W5 
complex may be a product of triggering, as it is located on 
one side of the W4 chimney thought to be created by mul- 
tiple supernovae during the last '-^lO MYr ( jOey et al.||2005 
Figure [T]) . 

Following [Koenig et al.] ([2008 ) , we adopt a distance to 
W5 of 2 kpc based on the water-maser parallax distance to 
the neighboring W3(0H) region ( [Hachisuka et al.|2006 ). As 
with W3, the W5 cloud is substantially 1.5 x) closer than 
its kinematic distance would suggest {vlsr{—^0 km s""*^) ~ 
3 kpc). Given this distance, Koenig et al. (2008) derived a 
total gas mass of 6.5x10^ Mq from a 2 fim extinction map. 

The W5 complex was mapped in the ^^CO 1-0 emis- 
sion line by the Five College Radio Astronomy Observatory 



(FCR AO) using the SEQUOIA receiver array ( Heyer et al. 
1998). The same array was used to map W5 in the ^"^CO 1- 



line (C. Brunt, private communication). Some early work 
searched for outflows in W5 ( [Bretherton et al.|2002 ), but the 
low-resolution CO 1-0 data only revealed a few, and only one 
was published. The higher resolution and sensitivity obser- 
vations presented here reveal many additional outflows. 

While W5 is thought to be associated with the W3/4/5 
complex, there are other infrared sources in the same part 
of the sky that are not obviously associated with W5. Some 
of these have been noted to be in the outer arm (several kpc 
behind W5) by[Digel et ar]([T996t and[Snell et"ar] ([2002t . 



In section 2, we present the new and archival data used 
in our study. In section 3, we discuss the outflow detection 
process and compare outflow detectability in W5 to that in 
Perseus. In section 4, we discuss the physical properties of 
the outflows and their implications for star formation in the 
W5 complex. In section 5, we briefly describe the outer- arm 
outflows discovered. 



2 OBSERVATIONS 

2.1 JCMT HARP CO 3-2 

CO J=3-2 345.79599 GHz data were acquired at the 15 m 
James Clerk Maxwell Telescope (JCMT) using the HARP 
array on a series of observing runs in 2008. On 2-4 January, 
2008, ^ 800 square arcminutes were mapped. During the 
run, T225, the zenith opacity at 225 GHz measured using 
the Caltech Submillimeter Observator tipping radiometer, 
ranged from 0.1 to 0.4 (0.4 < T345GHZ < l-^- Additional 
areas were mapped on 4-7 August, 16-20 and 31 October, 



^ [http : //docs . j ach . hawaii . edu/ JCMT/SCD/SN/002 . 2/node5 . 
|html| 



and 1 and 12-15 Nov 2008 in similar conditions. A total of ^ 
3 square degrees (12000 arcmin^) in the W5 complex were 
mapped (a velocity-integrated mosaic is shown in Figure |2|. 

HARP is a 16 pixel SIS receiver array acting as a front- 
end to the ACSIS digital auto-correlation spectrometer. In 
January 2008, 14 of the 16 detectors were functional. In 
the 2nd half of 2008, 12 of 16 were functional, necessitating 
longer scans to achieve similar S/N. 

In January 2008, a single spectral window centered at 
345.7959899 with bandwidth 1.0 GHz and channel width 488 
kHz (0.42 km s~"^) was used. In August 2008 and later, we 
used 250 MHz bandwidth and 61 kHz (0.05 km s~^) channel 
width. At this frequency, the beam FWHM is 14^^ (0.14 pc 
at a distance of 2 kpc) r\ 

A raster mapping strategy was used. In 2008, the array 
was shifted by 1/2 of an array spacing (58.2^^) between scans. 
Data was sampled at a rate of 0.6s per integration. Two 
perpendicular scans were used for each field observed. Most 
fields were 10x10^ and took ^ 45 minutes. When only 12 
receptors were available, 1/4 array stepping (29.1^^) was used 
with a sample rate of 0.4s per integration. 

Data were reduced using the SMURF package within 
the STARLINK software distribution The SMURF com- 
mand MAKECUBE was used to generate mosaics of contiguous 
sub-fields. The data were gridded on to cubes with 6" pix- 
els and smoothed with a cr = 2-pixel gaussian, resulting in 
a map FWHM resolution of 18^^(0.17 pc). A linear fit was 
subtracted from each spectrum over emission-free velocities 
(generally -60 to -50 and -20 to -10 km s"^) to remove the 
baseline. The final map RMS was ar^ ^ 0.06 — O.llK in 
0.42 km s~^ channels. 

The sky reference position (off position) in January 2008 
was J2000 2:31:04.069 +62:59:13.81. In later epochs, off po- 
sitions closer to the target fields were selected from blank 
sky regions identified in January 2008 in order to increase 
observing efficiency. A main-beam efficiency Tjmb = 0.60 was 
used as per the JCMT website to convert measurements to 
Tmb, though maps and spectra are presented in the original 
TX units. 



2.2 FCRAO Outer Galaxy Survey 

The FCRAO Out er Galaxy Survey (OGS) observed the W5 
complex in ^^CO ( |Heyer et al.|l998| ) and ^^CO 1-0 (C. Brunt, 
private communication). The CO data cube achieved a 
mean sensitivity of 0.35 K per 0.13 km s~^ channel, or 



0.6 K km s"' integrated. The '""CO cube was integrated 
over all velocities and resampled to match the BGPS map 
using the montage[^ package. The FWHM beam size was 
Ob =50'.'(0.48 pc). The integrated ^^CO data cube, with a 
sensitivity a — IK km s~^, is displayed with region name 
identifications in Figure |3] 



^ |http : / /docs . j ach . hawaii . edu/ JCMT/OVERVIEW/tell 
overview/ 

http : / / star link . j ach . hawaii . edu/ 
^ |http : //montage . ipac . caltech . edu/ 
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Figure 1. An overview of the W3/4/5 complex (also known as the "Heart and Soul Nebula") in false color. Orange shows 8 iim emission 
from the Spitzer and MSX satellites. Purple shows 21 cm continuum emission from the DRAO CGPS ( [Taylor et al.|20Q3| ; the DSS R 
image was used to set the display opacity of the 21 cm continuum as displayed (purely for aesthetic purposes). The green shows JCMT 
^■^CO 3-2 along with FCRAO -"^^CO 1-0 to fill in gaps that were not observed with the JCMT. The image spans ~ 7° in galactic longitude. 
This overview image shows the hypothesized interaction between the W4 superbubble and the W3 and W5 star-forming regions ( (Oey| 
|et al.|2005| >. 
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Figure 2. A mosaic of the CO 3-2 data cube integrated from -20 to -60 km s"-*^. The grayscale is hnear from to 150 K km s"-*^. The 
red and blue X's mark the locations of redshifted and blueshifted outflows. Dark red and dark blue plus symbols mark outflows at outer 
arm velocities. Green circles mark the location of all known BO and earlier stars in the W5 region from SIMBAD. 
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Figure 3. Individual region masks overlaid on the FCRAO ^^CO integrated image. The named regions, S201, AFGL4029, LWCas, 
W5NW, W5W, W5SE, W5S, and W5SW, were all selected based on the presence of outflows within the box. The inactive regions were 
selected from regions with substantial CO emission but without outflows. The 'empty' regions have essentially no CO emission within 
them and are used to place limits on the molecular gas within the east and west 'bubbles'. W5NWpc is compared directly to the Perseus 
molecular cloud in Section [3.1.11 



2.3 Spitzer 



Spitzer IRAC and MIPS 24 /xm images from |Koenig et aL] 
(2008) were used for morphological comparison. The reduc- 



tion and extraction techniques are detailed in their paper. 



3 ANALYSIS 

3.1 Outflow Detections 

Outflows were identified in the CO data cubes by manu- 
ally searching through position- velocity space for line wings 
using STARLINK's GAIA display software. Outflow candi- 
dates were identified by high velocity wings inconsistent with 
the local cloud velocity distribution, which ranged from a 
width of 3 km s~^ to 7 km s~^. Once an outflow candidate 
was identified in the posit ion- velocity diagrams, the velocity 
range over which the wing showed emission in the position- 
velocity diagram (down to TX = 0) was integrated over to 
create a map from which the approximate outflow size and 
position was determined (e.g. Figures [4land[5t 



Unlike [Curtis et al] ([2010) and Hatchell &: Dunham 
(2009), we did not use an 'objective' outflow identifica- 



tion method because of the greater velocity complexity and 
poorer spatial resolution of our observations. The outflow 
selection criteria in these papers requires the presence of 
a sub-mm clump in order to identify a candidate driving 
source (and therefore a targeted region in which to search for 
outflows), making a similar objective identiflcation impossi- 
ble for our survey. As discussed later in Section [4] the re- 
gions associated with outflows have wide lines and many are 
double-peaked. Additionally, many smaller areas associated 
with outflows have collections of gaussian-profiled clumps 
that are not connected to the cloud in position- velocity dia- 
grams but are not outflows. In particular, W5 is pockmarked 
by dozens of small cometary globules that are sometimes 
spatially coincident with the clouds but slightly offset in ve- 
locity. 



While Arce et al. (2010) described the beneflts of 3D 



visualization using isosurface contours, we found that the 



varying signal-to-noise across large-scale (^ 500 pixeP) re- 
gions with signiflcant extent in RA/Dec and limited velocity 
dynamic range made this method diffult for W5. There were 
many low-intensity outflows that were detectable by care- 
ful searches through position- velocity space that are not as 
apparent using isosurface methods. Out of the 55 outflows 
reported here, only 14 would be considered obvious, high- 
intensity, high-velocity flows from their spectra alone; the 
rest could not be unambiguously detected without a search 
through position- velocity space. 

In the majority of sources, the individual outflow lobes 
were unresolved, although some showed hints of position- 
velocity gradients at low significance and in many the red 
and blue fiows are spatially separated. Only Outfiow I's 
lobes were clearly resolved (Figure [4|. Some of most sug- 
gestive gradients occurred where the outfiow merged with 
its host molecular cloud in posit ion- velocity space, making 
the gradient difficult to distinguish (e.g., Outfiow 12, Figure 
|6|. Bipolar pairs were selected when there were red and blue 
fiows close to one another. The classification of a bipolar fiow 
was either 'yc' (yes - confident), 'yu' (yes - unconfident), or 
'n' (no) in Table [5] This identification is discussed in the 
captions for each outfiow figure in the online supplement. 
The AFGL 4029 region has many red and blue lobes but 
confusion prevented pairing. 

In cases where only the red- or blue-shifted lobe was 
visible, the surrounding pixels were searched for lower- 
significance and lower- velocity counterparts. For cases in 
which emission was detected, a candidate counterfiow was 
identified and incorporated into the catalog. However, in 12 
cases, the counterfiow still evaded detection, either because 
of confusion or because the counterfiow is not present in CO. 

The outfiow positions are overlaid on the CO 3-2 image 
in Figure |2] to provide an overview of where star formation 
is most active. The figures in Section [5 .2| show outfiow loca- 
tions overlaid on small-scale images. 



5 Outflows 15, 20, 24, the cluster of outflows 26-32, 47, 48, 52, 
and 53 could all have readily been detected by pointed single-dish 
measurements. 
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Figure 4. Position- velocity diagrams (a), spectra (c), and contour overlays of Outflow 1 on Spitzer 4.5 /im (b) and 8 iim (d) images. 
This outflow is clearly resolved and bipolar, (a): Position- velocity diagram of the blue flow displayed in arcsinh stretch from =0 to 
3 K. Locations of the red and blue flows are indicated by vertical dashed lines. The location of the position-velocity cut is indicated 
by the orange dashed line in panels (b) and (d), although the position- velocity cut is longer than those cut-out images, (b) Spitzer 4.5 
jj^m image displayed in logarithmic stretch from 30 to 500 MJy sr"-*^. (c): Spectrum of the outflow integrated over the outflow aperture 
and the velocity range specifled with shading. The velocity center (vertical dashed line) is determined by fltting a gaussian to the -"^^CO 
spectrum in an aperture including both outflow lobes. In the few cases in which 13 CO 1-0 was unavailable, a gaussian was flt to the 
■"^"^CO 3-2 spectrum, (d): Contours of the red and blue outflows superposed on the Spitzer 8 fim image displayed in logarithmic stretch. 
The contours are generated from a total intensity image integrated over the outflow velocities indicated in panel (c). The contours in 
both panels (b) and (d) are displayed at levels of 0.5,1,1.5,2,3,4,5,6 K km s"-*^ (cr ^ 0.25 K km s"-*^). The contour levels and stretches 
specifled in this caption apply to all of the flgures in the supplementary materials except where otherwise noted. 
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Figure 5. Position-velocity diagram, spectra, and contour overlays of Outflow 2 (see Figure[4]for a complete description). While the two 
lobes are widely separated, there are no nearby lobes that could lead to confusion, so we regard this pair as a reliable bipolar outflow 
identiflcation. 



Because our detection method involved searching for 
high-velocity outflows by eye, there should be no false de- 
tections. However, it is possible that some of these outflows 
are generated by mechanisms other than protostellar jets 
and winds since we have not identified their driving sources. 

One possible alternative driving mechanism is a photoe- 
vaporation flow, which could be accelerated up to the sound 
speed of the ionized medium, cu ^ 10 km s~^. Gas accel- 
erating away from the cloud would not be detected as an 
outflow because it would be rapidly ionized. However, gas 
driven inward would be accelerated and remain molecular. It 
could exhibit red and / or blue flows depending on the line 
of sight orientation. While there are viable candidates for 
this form of outflow impersonator, such flows can only have 
peak velocities v < c///4 ^ 2.5 km s~^ in the strong adia- 
batic shock limit, so that any gas seen with higher velocity 
tails are unlikely to be radiation-driven. 

Another plausible outflow impostor is the high- velocity 
tail in a turbulent distribution. However, for a typical molec- 
ular cloud, the low temperatures would require very high 
mach-number shocks (M > 10 assuming Tdoud ^ 20 K 
and Vfiow ^ 3 km s~^) that in idealized turbulence should 



be rare and short-lived. It is not known how frequent such 
high- velocity excursions will be in non-ideal turbulence with 
gravity (A. Goodman, P. Padoan, private communication). 
Finally, it is less likely for turbulent intermittency to have 
nearly coincident red and blue lobes, so intermittency can 
be morphologically excluded in most cases. 



3.1.1 Comparison to Perseus CO 3-2 observations 



We used the HARP CO 3-2 cubes from Hatchell et al. (2007) 



to evaluate our ability to identify outflows. We selected an 
outflow that was well-resolved and unconfused, LI 448, and 
evaluated it at both the native sensitivity of the |Hatchell| 
et al. (2007) observations and degraded in resolution and 



sensitivity to match our own. We focus on L1448 IRS2, la- 
beled Outflow 30 in [Hatchell et al.| ( |2007t . Figure |7| shows a 
comparison between the original quality and degraded data. 

Integrating over the outflow velocity range, we measure 
each lobe to be about 1.6' x 0.8' (0.14 x 0.07pc). Assum - 
ing a distance to Perseus of 250 pc (e.g. Enoch et al.|2006' ), 
we smooth by a factor of 8 by convolving the cube with 
a FWHM = 111'' gaussian, then downsample by the same 
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Figure 6. Position-velocity diagram, spectra, and contour overlays of Outflow 12. Much of the red outflow is lost in the complex velocity 
proflle of the molecular cloud(s), but it is high enough velocity to still be distinguished. 



factor of 8 to achieve 6'' square pixels at 2 kpc. The result- 
ing noise was reduced because of the spatial and spectral 
smoothing and was measured to be ~ 0.05 K in 0.54 km s~^ 
channels, which is comparable to the sensitivity in our sur- 
vey. It is still possible to distinguish the outflows from the 
cloud in position-velocity space. Each lobe is individually 
unresolved (long axis ^ V2!' compared to our beam FWHM 
of 18^^), but the two are separated by > 20^^ and therefore 
an overall spatial separation can still be measured. Because 
they are just barely unresolved at this distance, the lobes' 
surface brightnesses are approximately the same at 2 kpc as 
at 250 pc; if this outflow were seen at a greater distance it 
would appear fainter. 

Hatchell et al. ( 2007| detected 4 outflows within this 



map, plus an additional confused candidate. We note an ad- 
ditional grouping of outflowing material in the north-middle 
of the map(centered on coordinate 150x150 in Figure|7|. In 
the smoothed version, only three outflows are detected in the 
blue and two in the red, making flow- counter flow association 
difficult. The north-central blueshifted component appears 
to be the counterpart of the red flow when smoothed, al- 
though it is clearly the counterpart of the northwest blue 
flow in the full-resolution image. 



We are therefore able to detect any outflows compara- 
ble to L1448 (assuming a favorable geometry), but are likely 
to see clustered outflows as single or possibly extended lobes 
and will count fewer lobes than would be detected at higher 
resolution. Additionally, it is clear from this example that 
two adjacent outflows with opposite polarity are not neces- 
sarily associated, and therefore the outflows' source(s) may 
not be between the two lobes. 

In order to determine overall detectability of outflows 



compared to Perseus, we compare to [Curtis et al.| ( |2010| ) m 
Figure [8] Out of 29 outflows in their survey with measured 
'lobe lengths', 22 (71%) were smaller than 128^^ which would 
be below our 18^^ resolution if observed at 2 kpc. Even the 
largest lobes (HRF26R, HRF28R, HRF44B) would only ex- 
tend ^ 60^^ at 2 kpc. Each lobe in the largest outflow in our 
survey. Outflow 1, is ^ 80^^ (660^^ at 250pc), but no other in- 
dividual outflow lobes in W5 are clearly resolved. However, 
as seen in Figure [S] many bipolar lobes are separated by 
more than the telescope resolution, and the overall lobe sep- 
aration distribution (as opposed to the lobe length, which 
is mostly unmeasured in our sample) in W5 is quite similar 
to the separation distribution in Perseus. The 2-sample KS 
test gives a 25% probability that they are drawn from the 
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Figure 7. Comparison of L1448 seen at a distance of 250 pc (left) versus 2 kpc (middle) with sensitivity 0.5 K and 0.05K per 0.5 km 
channel respectively. Far Left Position-velocity diagram (log scale) of the outflow L1448 IRS2 at its native resolution and velocity. L1448 
IRS2 is the rightmost outflow in the contour plots. The PV diagram is rotated 45° from RA/Dec axes to go along the outflow axis. Middle 
Left: Posit ion- velocity diagram (log scale) of the same outflow smoothed and rebinned to be eight times more distant. Top Right: The 
integrated map is displayed at its native resolution (linear scale). The red contours are of the same data integrated from 6.5 to 16 km 
and the blue from -6 to km s"-*^. Contours are at 1,3, and 5 K km s"-*^ (~ 6, 18, 30cr). Axes are offsets in arcseconds. Because we are only 
examining the relative detect ability of outflows at two distances, we are not concerned with absolute coordinates. Bottom Right: The 
same map as it would be observed at eight times greater distance. Axes are offsets in arcseconds assuming the greater distance. Contours 
are integrated over the same velocity range as above, but are displayed at levels 0.25,0.50,0.75,1.00 K km (~ 12, 24, 48, 60cr). The 
entire region is detected at high signiflcance, but dominated by confusion. It is still evident that the red and blue lobes are distinct, but 
they are each unresolved. 
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Figure 8. Histogram of the measured outflow lobe separations. The grey hatched region shows [Curtis et al. [ ( 120101 values. The vertical 
dashed line represents the spatial resolution of our survey. The two distributions are similar. 
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same distribution (the null hypothesis that they are drawn 
from the same distribution cannot be rejected). 

On average, the Curtis et al. (2010) outflow velocities 
are similar to ours (Figure [9]). We detect lower velocity out- 
flows because we do not set a strict lower velocity limit cri- 
terion. We do not detect the highest velocity outflows most 
likely because of our poorer sensitivity to the faint high- 
velocity tips of outflows, although it is also possible that no 
high- velocity {v > 20 km s~^) flows exist in the W5 region. 
Note that the histogram compares quantities that are not 



directly equivalent: the outflows in Curtis et al. (2010) and 



our own data are measured out to the point at which the 
outflow signal is lost, while the 'region' velocities are full- 
width half-max (FWHM) velocities. 

Finally, we use the detectability of outflows in Perseus 
to inform our expectations in W5. Since it appears that 
we can detect outflows from low-mass protostars with sub- 
stellar to ^ 30Lq luminosities at the distance of W5 and 
these objects should be the most numerous in a standard 
initial mass function, the distribution of physical proper- 
ties in W5 outflows should be similar to those in Perseus. 
However, because W5 is a somewhat more massive cloud 
{Mw5 ~ ^Mperseus wc might cxpcct the high-end of 
the distribution to extend to higher values of outflow mass, 
momentum, and energy. Since we will likely see clustered 
outflows confused into a smaller number of distinct lobes, 
we expect a bias towards higher values of the derived quan- 
tities but a lower detection rate. 



3.1.2 Velocity, Column Density, and Mass Measurements 

Throughout this section, we assume that the CO lines are 
optically thin and thermally excited. The measured proper- 
ties are presented in Table [2] These assumptions are likely 
to be invalid, so we also discuss the consequences of applying 
'typical' optical depth corrections to the derived quantities. 
Because we do not measure optical depths and the optical 
depth correction for CO 3-2 is less well quantified than for 
CO 1-0 dCurtis et al.||2010| [Cabrit Bertout||l990f | we 
only present the uncorrected measurements in Table [3] 

The outflow velocity ranges were measured by exam- 
ining both RA-velocity and Dec-velocity diagrams interac- 
tively using the STARLINK GAIA data cube viewing tool. 
The velocity limits are set to include all outflow emission 
that is distinguishable from the cloud (i.e. the velocity at 
which outflow lobes dominate over the gaussian wing of the 
cloud emission) down to zero emission. An outflow size (or 
lobe size, following ^Curtis et al.||2010 ) was determined by 
integrating over the blue and red velocity ranges and creat- 
ing an elliptical aperture to include both peaks; the position 



^ M\Y5 is estimated from -"^"^CG. We also estimate the total molec- 
ular mass in W5 using the X-factor and acquire M\y5 = 5.0 x 10^ 
Mq, in agreement with |Karr &: Martin i (|2003t, who estimated a 
molecular mass of 4.4x10"^ from ^^CO using the same X-factor. 
[Koenig et al-l i l^OOSl estimated a total gas mass of 6.5x10^ from 
a 2MASS extincti on map. The tota l molecular mass in Perseus is 
Mpe rseus ^ 10^ ( Bally et al.|2008| > 

^ In [Curtis et al. (|2U1U|, this correction factor ranged from 1.8 



to 14.3; |Arce et al.| ( |2010| ) did not enumerate the optical depth 
correction they used but it is typically around 7 ( ^Cabrit Bertout, 
T990l >. 



and size therefore have approximately beam-sized (~ 18^^) 
accuracy. The integrated outflow maps are shown as red and 
blue contours in Figure |5] The velocity center was computed 
by fitting a gaussian to the FCRAO ^^CO spectrum aver- 
aged over the elliptical aperture. 

The column density is estimated from ^^CO J=3- 
2 assuming local thermal equilibrium (LTE) and opti- 
cally thin emission using the equation A/'(H2) = 5.3 x 
10^^r]~l f TX{v)dv for Tex = 20 K. The derivation IS given in 
the Appendix. The column density in the lobes is likely to be 
dominated by low-velocity gas and therefore our dominant 
uncertainty may be missing low-velocity emission rather 
than poor assumptions about the optical depth. 

The scalar momentum and energy were computed from 



p = M 



J:TX{v){v - Vc)Av 
J:tx(v)Av 



ETXiv)Av 



(1) 



(2) 



where Vc is the ^"^CO 1-0 centroid velocity. The same as- 
sumptions used in determining column density are applied 
here. 

We estimate an outflow lifetime by taking half the dis- 
tance between the red and blue outflow centroids divided 
by the maximum measured velocity difference {Avmax = 

{jJraax,red max , blue) / •> "Tflow — L flow / i^^'^raax^ •> whcrC 

Lflow refers to the length of the flow. This method assumes 
that the outflow inclination is 45°; if it is more parallel to 
the plane of the sky, we overestimate the age, and vice- versa. 
The momentum flux is then P — p/r. Similarly, we compute 
a mass loss rate by dividing the total outflow mass by the 
dynamical age, which yields what is likely a lower limit on 
the mass loss rate (if the lifetime is underestimated, the mass 
loss rate is overestimated, but the outflow mass is always a 
lower limit because of optical depth and confusion effects). 

The dynamical ages are highly suspect since the red and 
blue lobes are often unresolved or barely resolved, and dif- 
fuse emission averaged with the lobe emission can shift the 
centroid position. Additionally, it is not clear what portion 
of the outflow corresponds to the centroid: the bow shock 
or the jet could both potentially dominate the outflow emis- 
sion. Curtis et al. (2010) discuss the many ways in which 



the dynamical age can be in error. Our mass loss rates are 
similar to those in Perseus without correcting our measure- 
ments for optical depth, while our outflow masses are an 
order of magnitude lower. It therefore appears that our dy- 
namical age estimates must be too low, since we have no 
reason to expect protostars in W5 to be undergoing mass 
loss at a greater rate than those in Perseus. However, given 
more reliable dynamical age estimates from higher resolu- 
tion observations of shock tracers, the mass loss rates could 
be corrected and compared to other star-forming regions. 

Because the emission was assumed to be optically thin, 
the mass, column, energy, and momentum measurements we 
present are strictly lower limits. While some authors have 
computed correction factors to ^^CO 1-0 optical depths (e.g. 
Cabrit Bertout|1990 ), the corrections are different for the 
3-2 transition (1.8 to 14.3, [Curtis et al.|2010| ). Additionally, 
CO 3-2 may require a correction for sub-thermal excitation 
because of its higher critical density (the CO 3-2 critical 
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Figure 9. Histogram of the outflow line widths. Black lines: histogram of the measured outflow widths (half- width zero-intensity, 
measured from the fltted central velocity of the cloud to the highest velocity with non-zero emission). Blue dashed lines: outflow half- 
width zero-intensity (HWZI) for the outer arm (non-W5) sample. Solid red shaded: The measured widths (HWHM) of the sub-regions 
as tabulated in Table [T| Gray dotted: Outflow Vmax 

(HWZI) values for Perseus from |Curtis et al.|p010t . 



density is 27 times higher than CO 1-0; see Appendix [A] for 
modeling of this effect). 

Additionally, most of the outflow mass is at the lowest 
distinguishable velocities in typical outflows (e.g.jArce et al. 



2010 ). It is therefore plausible that in the more turbulent W5 



region, a greater fraction of the outflow mass is blended (ve- 
locity confused) with the cloud and therefore not included 
in mass, momentum, and energy measurements. This omis- 
sion could be greater than the underestimate due to poor 
opacity assumptions. 

The total mass of the W5 outflows is Mtot ~ l.SM©, 
substantially lower, even with an optical depth correction 
of 10 X, than the 163 Mq reported in Perseus (Arce et al. 



2010). Arce et al. (2010) also include a correction factor of 



2.5 to account for higher temperatures in outflows and a fac- 
tor of 2 to account for emission blended with the cloud. The 
temperature correction is inappropriate for CO 3-2 (see Ap- 
pendix [A] Figure Al), but the resulting total outflow mass 
in W5 with an optical depth correction and a factor of 2 
confusion correction is about 30 Mq. In order to make our 
measurements consistent with a mass of 160 Mq , a density 
upper limit in the outflowing gas of n(H2) < 10^"^cm~^ is 
required, since a lower gas density results in greater mass 
for a given intensity (see Appendix [A] Figure A2 ). However, 
we expect the total outflow mass in W5 to be greater than 
in Perseus because of the greater cloud mass, implying that 
the density in the flows must be even lower, or additional 
corrections are needed. 

The total outflow momentum is ptot ^ 10.9Mq km s~^, 



versus a quoted 517 Mokm s ^ in Perseus (Arce et al. 2010 ) 



Arce et al. ( 2010p included inclination and dissociative shock 
corrections for the momentum measurements on top of the 
correction factors already applied to the mass. If these cor- 
rections are removed from the Perseus momentum total (ex- 
cept for optical depth, which is variable in their data and 
therefore cannot be removed), the uncorrected outflow mo- 
mentum in Perseus would be about 74 Mokm s~^. The 
W5 outflow momentum, if corrected with a 'typical' opti- 
cal depth in the range 7-14, would match or exceed this 



value. If an additional CO 3-2 excitation correction (in the 
range 1-20) is applied, the W5 outflow momentum would 
significantly exceed that in Perseus. 

Assuming a turbulent line width Av ^ 3 km s~^ (ap- 
proximately the smallest FWHM line- width observed), the 
total turbulent momentum in the ambient cloud is p — 
MtotAv = 1.3 X lO^M© km s~\ which is 10^ times the 
measured outflow momentum - the outflows detected in our 
survey cannot be the sole source of the observed turbulent 
line widths, even if corrected for optical depth and missing 
mass. 

Table ^presents the turbulent momentum for each sub- 
region computed by multiplying the measured velocity width 
by the integrated ^^CO mass. Even if the outflow measure- 
ments are orders of magnitude low because of optical depth, 
cloud blending, sub-thermal excitation, and other missing- 
mass considerations, outflows contribute negligibly to the 
total momentum of high velocity gas in W5. This result is 
unsurprising, as there are many other likely sources of en- 
ergy in the region such as stellar wind bubbles and shock 
fronts between the ionized and molecular gas. Additionally, 
in regions unaffected by feedback from the HII region (e.g. 
W5NW), cloud-cloud collisions are a possible source of en- 
ergy. 

Figure [To] displays the distribution of measured proper- 
ties and compares them to those derived in the COMPLETE 

( [20T0 t HARP CO 3- 



(Arce et al. 2010) and Curtis et al 



2 surveys of Perseus. Our derived masses are substantially 
lower than those in Arce et al. (2010) even if corrected for 



optical depth, but our momenta are similar to the CPOC 
(COMPLETE Perseus Outflow Candidate) sample and our 
energies are higher, indicating a bias towards detecting mass 
at high velocities. The bias is more heavily towards high ve- 
locities than the CO 1-0 used in "Arce et al.|(|2010|). The dis- 



crepancy between our values and those of Arce et al. (2010) 



and Curtis et al. (2010) can be partly accounted for by the 



optical depth correction applied in those works: ^^CO was 
used to correct for opacity at low velocities, where most of 
the outflow mass is expected. Those works may also have 
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Figure 10. Histograms of outflow physical properties. The sohd unfllled hnes are the W5 outflows (this paper), the forward-slash hashed 
lines show [Arce et al.| ( |2010[ ) CPOCs , the dark gray shaded region shows [Arce et al!] ( |2010| ) values for known outflows in Perseus, and 
the light gray, backslash-hashed regions show |Curtis et ah] ( |2010| CO 3-2 outflow properties. The outflow masses measured in Perseus 
are systematically higher partly because both surveys corrected for line optical depth using -"^^CO. The medians of the distributions are 
0.017, 0.044, 0.33, and 0.14 Mq for W5, Curtis, Arce Known, and Arce CPOCs respectively, which implies that an optical depth and 
excitation correction factor of 2.5-20 would be required to make the distributions agree (although W5, being a more massive region, 
might be expected to have more massive and powerful outflows). It is likely that CO 3-2 is sub-thermally excited in outflows, and CO 
outflows may be destroyed by UV radiation in the W5 complex while they easily survive in the lower-mass Perseus region, which are 
other factors that could push the W5 mass distribution lower. 



been less affected by blending because of the smaller cloud 
line widths in Perseus. 

The momentum flux and mass loss rate are compared 



to the values derived in Perseus by "Hatchell et aL| (|2007 ) 
and |Curtis et aL| (|2010| in Figures 



11 



and 



12 



Both of our 



values are computed using the dynamical timescale Td mea- 
sured from outflow lobe separation, while the [Hatchell et dl\ 
(2007) values are derived using a more direct momentum- 
flux measurement in which the momentum flux contribution 
of each pixel in the resolved outflow map is considered. The 



derived momentum fluxes (Figure 11) are approximately 
consistent with the|Curtis et al. (2010) Perseus momentum 



fluxes; |Curtis et al 



(2010) measure momentum fluxes in a 
<P < 7x 10"^ Mokms 



yr 



higher than 



range 1 x 10 

our measured 6 x 10~^ < P < 1 x 10~^ Mokm s ^yr 
by approximately the opacity correction they applied. As 
seen in Figure |11| the Hatchell et al. (2007) momentum 
flux measurements in Perseus cover a much lower range 
6 X 10"^ < P < 2 X 10"^ Mokm s~^yr~^ and are not con- 



sistent with our measurements. This disagreement is most 
likely because of the difference in method. The W5 outflows 
exhibit substantially higher mass-loss rates and momentum 
fluxes if we assume a factor of 10 opacity correction, as 
expected from our bias toward higher- velocity, higher-mass 
flows. 



3.2 Structure of the W5 molecular clouds: A thin 
sheet? 

The W5 complex extends - 1.6° x 0.7° within 20° of par- 
allel with the galactic plane. At the assumed distance of 
2 kpc, it has a projected length of ^ 60 pc (Figure |2|. 
In the 8 lira band (Figure [T]), the region appears to con- 
sist of two blown-out bubbles with ^ 10 — 15 pc radii cen- 
tered on ^ = 138.1,6 = 1.4 and ^ = 137.5,6 = 0.9. While 
the bubbles are filled in with low-level far-infrared emission, 
there is no CO detected down to a 3 — a limit of 3.0 K 
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Figure 11. Histogram of the measured outflow momentum fluxes. The black thick hue shows our data, the grey shaded region shows the 
[Hatchell et al.| ( |2007| > data, and the hatched region shows Curt is et al.| ( |2010| ) values. Our measurements peak squarely between the two 
Perseus JCMT CO 3-2 data sets, although the Curtis et al. (2Q10| results include an opacity correction that our data do not, suggesting 
that our results are likely consistent with |Curtis et al.| ( |2010| ) but inconsistent with the [Hatchell et al.| ( |2007j ) direct measurement method. 




Figure 12. H istogram of the measured mass loss rate. The black thick line shows our data, while the grey shaded region shows the 
Hatchell et al. (2007) data, which is simply computed by M = P x 10/5 km s"-*^, where the factor of 10 is a correction for opacity. Our 
mass loss rates are very comparable to those of Hat chell et al.| ( |2007| , but different methods were used so the comparison may not be 
physically meaningful. [Curtis et al.| ( |2Q1Q| ) (hatched) used a dynamical time method similar to our own and also derived similar mass 
loss rates, although their mass measurements have been opacity-corrected using the ^^CO 3-2 line. Because our mass loss rates agree 
reasonably with Perseus, but our outflow mass measurements are an order of magnitude low, we believe our dynamical age estimates to 
be too small. 
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km s-^ (^^CO 1-0), 2.4 K km s"^ (^^CO 3-2, excepting a 
few isolated clumps), and 1.5 K km s"^ (^^CO 1-0). Us- 
ing the X-factor (the CO-to-H2 conversion factor) for ^^CO 
iV(H2) = 3.6 X 10^°cm~^/(K km s~^), we derive an upper 
limit A^(H2) < 1.1 x lO^^cm"^ or Ay < 0.6. Individual 
'wisps' and 'clumps' of CO can sometimes be seen, partic- 
ularly towards the cloud edges, but in general the bubbles 
are absent of CO gas. 

Given such low column limits, the W5 cloud must be 
much smaller along the line of sight than its ~ 50 pc size 
projected on the sky. Alternately, along the line-of-sight, 
the columns of molecular gas are too low for CO to self- 
shield, and it is therefore destroyed by the UV radiation 
of W5's 0-stars. In either case, there is a significant excess 
of molecular gas in the plane of the sky compared to the 
line of sight, which makes W5 an excellent location to per- 
form unobscured observations of the star formation process. 
The implied thin geometry of the W5 molecular cloud may 
therefore be similar to the bubbles observed by IBeaumo'nt] 
& Williams (2010), but on a larger scale. 



There is also morphological evidence supporting the 
face-on hypothesis. In the AFGL 4029 region (Section 4.2) 
and all along the south of W5, there are ridges with many 
individual cometary 'heads' pointing towards the 0-stars 
that are unconfused along the line of sight. This sort of sep- 
aration would not be expected if we were looking through 
the clouds towards the 0-stars. W5W, however, presents a 
counterexample in which there are two clouds along the line 
of sight that may well be masking a more complex geometry. 



4 SUB-REGIONS 

Individual regions were selected from the mosaic for com- 
parison. All regions with multiple outflows and indicators of 
star formation activity were named and included as regions 
for analysis. Additionally, three "inactive" regions were se- 
lected based on the presence of ^^CO emission but the lack 
of outflows in the ^^CO 3-2 data. Finally, two regions devoid 
of CO emission were selected as a baseline comparison and 
to place upper limits on the molecular gas content of the 
east and west 'bubbles'. The regions are identified on the 
integrated ^^CO image in Figure [s] 

Average spectra were taken of each "region" within the 
indicated box. Gaussians were fit to the spectrum to deter- 
mine line-widths and centers (Figure 13 Table [TJ. Gaussian 
fits were necessary because in many locations there are at 
least two velocity components, so the second moment (the 
"intensity- weighted dispersion" ) is a poor estimator of line 
width. Widths ranged from vfwhm — 2.3 to 6.2 km s~^ 
(Figure [9]). 



4.1 Sh 2-201 

Sh 2-201 is an HII region and is part of the same molecular 
cloud as the bright-rimmed clouds in W5E, but it does not 



share a cometary shape with these clouds (Figure 14). In- 
stead, it is internally heated and has its own ionizing source 
( Fein et al.|p^87] ). The AFGL 4029 cloud edge is at a pro- 
jected distance of ~ 7 pc from the nearest exposed 0-star, 
and the closest illuminated point in the Spitzer 8 and 24 iim 
maps is at a projected distance of 5 pc. The star forming 



process must therefore have begun before radiation driven 
shocks from the W5 0-stars could have impacted the cloud. 



4.2 AFGL 4029 

AFGL 4029 is a young cluster embedded in a cometary cloud 



(Figure 15). There is one clear bipolar outflow and 6 single- 
lobed flows that cannot be unambiguously associated with 
an opposite direction counterpart. The cluster is mostly un- 
resolved in the data presented here and is clearly the most 
active CO clump in W5. It contains a cluster of at least 30 
B-stars ( Deharveng et al.]|1997 ). The outflows from this re- 
gion have a full width A^' ^ 30 km s~^, which is entirely 
inconsistent with a radiation-driven inflow or outflow since 
it is greater than the sound speed in the ionized medium. 

The northeast cometary cloud is strongly affected by the 
W5 HII region. It has an outflow in the head of the cloud 
(Figure 16), and the cloud shows a velocity gradient with 



distance from the HII region. The polarity of the gradient 
suggests that the cometary cloud must be on the far side of 
the ionizing 0-star along the line of sight assuming that the 
HII region pressure is responsible for accelerating the cloud 
edge. 



4.3 W5 Ridge 

The W5 complex consists of two HII region bubbles sepa- 
rated by a ridge of molecular gas (Figure 17). This ridge con- 



tains the LW Gas optical nebula, a reflection nebula around 
the variable star LW Gas, on its east side and an X-shaped 
nebula on the west. The east portion of LW Gas Nebula is 
bright in both the continuum and CO J=3-2 but lacks out- 



flows (see Figure 17). The east portion also has the highest 
average peak antenna temperature, suggesting that the gas 
temperature in this region is substantially higher than in the 
majority of the W5 complex (higher spatial densities could 
also increase the observed Ta, but the presence of nearby 
heating sources make a higher temperature more plausible). 
It is possible that high gas temperatures are suppressing 
star formation in the cloud. Alternately, the radiation that 
is heating the gas may destroy any outflowing CO, which 
is more likely assuming the two Class I objects identified in 
this region by Koe nig et al.| (2008) are genuine protostars. 
The ridge is surprisingly faint in HI 21 cm emission 



compared to the two HII regions (Figure 18) considering its 
24 iim surface brightness. The integrated HI intensity from 
-45 to -35 km s~^ is ~ 800 K km s~\ whereas in the HII 
region bubble it is around 1000 K km s"\ The CO-bright 
regions show lower levels of emission similar to the ridge 
at 700-800 K km s~^. However, the ridge contains no CO 



gas and very few young stars (Figure 7 in [Koenig et al. 
2008 ) . It is possible that the ridge contains cool HI but has 



very low column-densities along the direction pointing to- 
wards the 0-stars, in which case the self-shielding is too lit- 
tle to prevent CO dissociation. This ridge may therefore be 
an excellent location to explore the transition from molec- 
ular to atomic gas under the influence of ionizing radiation 
in conditions different from high-density photodissociation 
(photon-dominated) regions. 

We examine Outflow 20 as a possible case for pressure- 
driven implosion (radiation, RDI, or gas pressure, PDI) 
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Table 1. Gaussian fit parameters of sub-regions 



Region 


Velocity 1 


VV IQIjII 1 


Amplitude 1 


veiociijy z 


VV IQIjII Z 


/\iiipiiijUQe Z 




(km s-l) 


(FWHM, km s-i) 


(K) 


(km s-l) 


(FWHM, km s-i) 


(K) 


S201 


-38.04 


3.149 


2.35 








AFGL4029 


-38.91 


3.3605 


1.48 








LWCas 


-38.83 


3.478 


2.33 








W5W 


-41.37 


3.8775 


3.07 


-36.16 


3.8305 


1.90 


W5NW 


-36.37 


3.854 


1.6 








WSNWpc 


-36.37 


3.713 


1.19 


-41.81 


4.3475 


0.47 


W5SW 


-42.78 


4.136 


0.6 


-36.34 


4.183 


0.22 


W5S 


-40.15 


2.914 


0.34 


-35.76 


2.2795 


0.40 


Inactive 1 


-42.91 


2.6555 


0.75 


-39.38 


4.2065 


0.42 


Inactive2 


-38.94 


3.7365 


1.2 








empty 


-37.81 


5.217 


0.04 









l^CO fits 13QO 

mass momentum energy 
(Mo) (Mokms-i) (ergs) 
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Figure 13. Spatially averaged spectra of the individual regions analyzed. -"^^CO 3-2 is shown by thick black lines and ^^CO 1-0 is shown 
by thin red lines. Gaussian fits are overplotted in blue and green dashed lines, respectively. The fit properties are given in Table [l] 
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Figure 14. Small scale map of the Sh 2-201 region plotted with CO 3-2 contours integrated from -60 to -20 km s"-*^ at levels 3,7.2,17.3,41.6, 
and 100 K km s"-*^. The IRAC 8 iim image is displayed in inverted log scale from 800 to 8000 MJy sr"-*^. Contours of the CO 3-2 cube 
integrated from -60 to -20 km s~ are overlaid at logarithmically spaced levels from 3 to 10 K km s-^ (3.0,7.2,17.3,41.6,100; cr 0.7 K 
km s~^). The ellipses represent the individual outflow lobe apertures mentioned in Section|3.1.2| 





; 1 ' , 
















JJIX- ■■".1 








1 


1 



I . \ 

133. a- Ti'^.A- 133.1- 

(jo aciic Long tudc 



Figure 15. Small scale map of the AFGL 4029 region plotted with CO 3-2 contours integrated from -60 to -20 km s -"^ at levels 
3,7.2,17.3,41.6, and 100 K km s-^ The IRAC 8 /im image is displayed in inverted log scale from 800 to 8000 MJy sr-^ Contours of the 
CO 3-2 cube integrated from -60 to -20 km are overlaid at logarithmically spaced levels from 3 to 100 K km s"-*^ (3.0,7.2,17.3,41.6,100; 
cr 0.7 K km s"-*^). Outflows 26-32 are ejected from a forming dense cluster. A diagram displaying the kinematics of the northern cometary 
cloud is shown in Figure 
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Figure 16. The northeast cometary cloud. Contours are shown at 0.5,1,2, and 5 K km s"-*^ integrated over the ranges -44.0 to -41.9 
km s"-*^ (blue) and -38.1 to -35.6 km (red). There is a velocity gradient across the tail, suggesting that the front edge is being pushed 
away along the line of sight. 
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Figure 17. Small scale map of the LW Cas nebula plotted with CO 3-2 contours integrated from -60 to -20 km s ~-^ at levels 3,7.2,17.3,41.6, 
and 100 K km s"-*^. The feature containing outflows 20 and 21 is the X-shaped ridge referenced in Section [4.3| This sub-region is notable 
for having very few outflows associated with the most signiflcant patches of CO emission. The gas around it is heated on the left side 
by the 07V star HD 18326 (Dproj = 8.5 pc), suggesting that this gas could be substantially warmer than the other molecular clouds in 
W5. 



by examining the relative timescales of the outflow driv- 
ing source and the Hll-region-driven compression front. A 
typical molecular out flow source (Class or I) has a lifetime 
of 5 X 10^ years (Evans et al. 



2009t . 



Given that there 

is an active outflow at the head of this cloud, we use 0.5 
MYr as an upper limit. The approximate distance from this 
source to the cloud front behind it is ^ 3.3 pc. If we as- 
sume the cloud front has been pushed at a constant speed 
V ^ cii ^ 10km s~^, we derive a lower limit on its age of 0.3 
MYr. While these limits allow for the protostar to be older 
than the compression front by up to 0.2 MYr, it is likely that 
the compression front moved more slowly (e.g., 3 km s~^ if 
it was pushed by a D-type shock front) and that the proto- 
star is not yet at the end of its lifetime - it is very plausible 
that this soure was born in a radiation-driven implosion. 



4.4 Southern Pillars 

There are 3 cometary clouds that resemble the "elephant 



contains evidence of at least one outflow in the head of the 
cloud (see the supplementary materials, outflows 16-19 and 
38) These pillars are low-mass and isolated; there is no other 
outflow activity in southern W5. However, because of the 
bright illumination on their northern edges and robust star 
formation tracers, these objects present a reasonable case 
for triggered star formation by the RDI mechanism. 

The kinematics of these cometary clouds suggest that 
they have been pushed in different directions by the HIT 
region (Figure 19). The central cometary cloud (Figure 19 3) 
has two tails. The southwest tail emission peaks around - 
39.5 km s~^ and the southeast tail peaks at -41.5 km s~^, 
while the head is peaked at an intermediate -40.5 km s~^. 
These velocity shifts suggest that the gas is being accelerated 
perpendicular to the head-tail axis and that the southeast 
tail is on the back side of the cometary head, while the 
southwest tail is on the front side. The expanding HIT region 
is crushing this head-tail system. 



trunk" nebula in IC 1396 (Figure 19). Each of these pillars 



The southeast cometary cloud (Figure 19 1) peaks at - 
35.0 km s"^ There are no clearly-separated CO tails as in 
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Figure 18. Top: The DRAO 21 cm HI map integrated from -45 to -35 km s"-*^ displayed in grayscale from 700 (black) to 1050 (white) 
K km s"-*^ with IRAS 100 /im contours (red, 40 MJy sr"-*^) and -"^^CO 1-0 contours integrated over the same range (white, 4 K km s"-*^) 
overlaid. The ridge of IRAS 100 fim emission at i = 138.0 coincides with a relative lack of HI emission at these velocities, suggesting 
either that there is less or colder gas along the ridge. Bottom: The Spitzer 24 /im map with 21 cm continuum contours at 6, 8, and 10 
MJy sr"-*^ overlaid. The IRAS contours are also overlaid to provide a reference for comparing the two figures and to demonstrate that 
the HII region abuts the cold-HI area. The moderate excess of continuum emission implies a somewhat higher electron density along the 
line of sight through the ridge. 
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the central cloud, but there is a velocity shift across the tail, 
in which the west (right) side is blueshifted compared to the 
east (left) side, which is the opposite sense from the central 
comet ary cloud. 

The southwest cometary cloud (Figure 19:) peaks at 



-40.3 km s~ and has weakly defined tails similar to the 
central cloud. Both of its tails are at approximately the same 
velocity (-42.5 km s~^). 

The kinematics of these tails provide some hints of their 
3D structure and location in the cloud. Future study to com- 
pare the many cometary flows in W5 to physical models and 
simulations is warranted. Since these flows are likely at dif- 
ferent locations along the line of sight (as required for their 
different velocities), analysis of their ionized edges may al- 
low for more precise determination of the full 3D structure 
of the clouds relative to their ionizing sources. 



4.5 W5 Southeast 

The region identified as W5SE has very little star for- 
mation activity despite having significant molecular gas 
(Misco 8OOM0). While there are two outflows and two 
Class I objects ( Koenig et al.|[2008 ) in the southeast of the 
two clumps {£ = 138.15,5 = 0.77, Figure [20|, the main 
clump {£ = 138.0,6 = 0.8) has no detected outflows. The 
CO emission is particularly clumpy in this region, with many 
independent, unresolved clumps both in position and veloc- 
ity. In the 8 and 24 micron Spitzer images, it is clear that 
these clouds are illuminated from the northwest. This region 
represents a case in which the expanding HII region has im- 
pacted molecular gas but has not triggered additional star 
formation. The high clump-to-clump velocity dispersion ob- 
served in this region may be analogous to the W5S cometary 



clouds (Section 4.4) but without condensed clumps around 



which to form cometary clouds. 



4.6 W5 Southwest 

There is an isolated clump associated with outflows in the 
southwest part of W5 (Figure 21) at vlsr —45 km s~"^. 
While this clump is likely to be associated with the W5 
region, it shows little evidence of interaction with the HII 
region. If it is eventually impacted by the expanding ioniza- 
tion front (i.e. if it is within the W5 complex), this clump 
will be an example of "revealed" , not triggered, star forma- 
tion. 

The other source in W5SW is a cometary cloud with 



a blueshifted head and redshifted tail (Figure 22 Outflow 
13). The head contains a redshifted outflow; no blueshifted 
counterpart was detected (the velocity gradient displayed in 
Figure [22] is smaller than the outflow velocity and is also 
redshifted away from the head). The lack of a blueshifted 
counterpart may be because the flow is blowing into ionized 
gas where the CO is dissociated. 

Because of its evident interaction with the HII re- 
gion, this source is an interesting candidate for a non- 
protostellar outflow impersonator. However, because the 
head is blueshifted relative to the tail, we can infer that 
the head has been accelerated towards us by pressure from 
the HII region, implying that it is in the foreground of the 
cloud. Given this geometry, a radiation-driven flow would 



appear blueshifted, not redshifted, as the detected flow is. 
Additionally, the outflow is seen as fast as 7.5 km s~^ red- 
shifted from the cloud, which is a factor of 2 too fast to be 
driven by radiation in a standard D-type shock. Finally, the 
outflow velocity is much greater than seen in a simulation 



of a cometary cloud by Gritschneder et al. (2010), while the 



head- to- tail velocity gradient is comparable. 



4.7 W5 West / IC 1848 

There is a bright infrared source seen in the center of W5W 
(IRAS 02459+6029; Figure |23]), but the nearest CO outflow 
lobe is 1 pc away. The nondetection may be due to con- 
fusion in this area: there are two layers of CO gas separated 
by ~5 km s~^, so low- velocity outflow detection is more 
difficult. Unlike the rest of the W5 complex, this region ap- 
pears to have multiple independent confusing components 



along the line of sight (Figure 13), and therefore the CO 



data provide much less useful physical information (multi- 
ple components are also observed in the ^"^CO data, ruling 
out self-absorption as the cause of the multiple components). 



4.8 W5 NW 

The northwest cluster containing outflows 1-8 is at a slightly 
different velocity (^ —35 km s~^) than the majority of the 
W5 cloud complex (^ —38 km s~^; Figure 24), but it shares 



contiguous emission with the neighboring W5W region. It 
contains many outflows and therefore is actively forming 



stars (Figure 25). However, this cluster exhibits much lower 
CO brightness temperatures and weaker Spitzer 8 /jm emis- 
sion than the "bright-rimmed clouds" seen near the W5 O- 
stars. We therefore conclude that the region has not been 
directly impacted by any photoionizing radiation from the 
W5 0-stars. 

The lack of interaction with the W5 0-stars implies 
that the star formation in this region, though quite vigor- 
ous, has not been directly triggered. Therefore not all of the 
current generation of star formation in W5 has been trig- 
gered on small or intermediate scales (e.g., radiation-driven 
implosion). Even the "collect and collapse" scenario seems 
unlikely here, as the region with the most outflows also dis- 
plays some of the smoothest morphology (Figures [2] and 25); 
in "collect and collapse" the expansion of an HII region leads 
to clumping and fragmentation, and the spaces between the 
clumps should be relatively cleared out. 



5 DISCUSSION 

5.1 Comparison to other outflows 

The outflow properties we derive are similar to those in the 
BO-star forming clum p IRAS 05358+3543 (M ^ GOOM© 
Ginsburg et al.|2009 ), in which CO 3-2 and 2-1 were used to 
derive outflow masses in the range 0.01-0.09 Mq. However, 
some significantly larger outflows, up to 1.6 pc in one direc- 
tion were detected, while the largest resolved outflow in our 
survey was only 0.8 pc (one direction). 

As noted in Section |3.1.1| the total molecular mass 
in W5 is larger than Perseus, Mwb 4.5 x IO^Mq while 
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Figure 19. CO 3-2 contours overlaid on the Spitzer 8 iim image of the W5S cometary clouds described in Section [4!4] Contours are 
color-coded by velocity and shown for 0.84 km channels at levels of 1 K (a, b) and 0.5 K (c). The velocity ranges plotted are (a) -41.5 
to -33.0 km s-i(b) -44.7 to -36.7 km s-^ (c) -43.6 to -35.6 km s-^ The labels show the minimum, maxmimum, and middle velocities to 
guide the eye. The grey boxes indicate the regions selected for CO contours; while there is CO emission associated with the southern 8 
lim emission, we do not display it here. The velocity gradients are discussed in Section |4.4| 
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Figure 20. Small scale map of the W5 SE region showing the star- forming clump containing outflows 39 and 40 and the non-star-forming 
clump at = 138.0, h = 0.8. CO 3-2 contours integrated from -60 to -20 km s-^ are displayed at levels 3,7.2,17.3,41.6, and 100 K km s-^ 
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Figure 21. Small scale map of the W5 SW region plotted with CO 3-2 contours integrated from -60 to -20 km s"-*^ at levels 3,7.2,17.3,41.6, 
and 100 K km s"-*^. Outflow 13 is at the head of a cometary cloud (Figure [22| and therefore has clearly been affected by the expanding 
HII region, but the region including bipolar Outflow 10 shows no evidence of interaction with the HIT region. 
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Figure 22. The cometary cloud in the W5 Southwest region (Outflow 13). Contours are shown at 1 K for 0.84 km s"-*^ wide channels 
from -37.2 km (blue) to -30.5 km s"-*^ (red). The head is clearly blueshifted relative to the tail and contains a spatially unresolved 
redshifted outflow. 
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Figure 23. Small scale map of the W5 W region. The IRAC 8 iim image is displayed in inverted log scale from 800 to 8000 MJy sr"-*^. 
Contours of the CO 3-2 cube integrated from -50 to -38 km s"-*^ (blue) and -38 to -26 km s"-*^ (red) are overlaid at levels 5,10,20,30,40,50,60 
K km a ^ 0.5 K km s"-*^. The lack of outflow detections is partly explained by the two spatially overlapping clouds that are adjacent 
in velocity. 
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Figure 24. Integrated longitude- velocity diagram of the W5 complex from b = 0.25 to 6 = 2.15 in -"^^CO 1-0 from the FCRAO OGS. 
The W5NW region is seen at a distinct average velocity around i = 136.5, v^sR = ~34 km s"-*^. The red and blue triangles mark the 
longitude-velocity locations of the detected outflows. In all cases, they mark the low-velocity start of the outflow. 
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Figure 25. Small scale map of the W5 NW region plotted with CO 3-2 contours integrated from -60 to -20 km s"-*^ at levels 3,7.2,17.3,41.6, 
and 100 K km s"-*^. Despite its distance from the W5 O-stars, Dproj ~ 20 pc, this cluster is the most active site of star formation in the 
complex as measured by outflow activity. 



M Perseus ^ 10^ M© ( [Bally et al.|2008| ). The length distribu- 
tion of outflows (Figure |8| is strikingly similar, while other 
physical properties have substantially different mean values 
with or without correction factors included. 

The W5NW region is more directly comparable to 
Perseus, with a total mass of - 1.5 X 10^ Mo (Table [l]) 
and a similar size. In Figure [3] we show both the W5NW 
region, which contains all of the identified outflows, and 
the WSNWpc region, which is a larger area intended to 
be directly comparable in both mass and spatial scale to 
the Perseus molecular cloud. The W5NWpc region contains 
more than an order of magnitude more turbulent energy 
than the Perseus complex {Eturb,Per = 1.6 x 10^^ ergs, 
Arce et aL||2010 ) despite its similar mass. Even the smaller 
W5NW region has 5x more turbulent energy than the 
Perseus complex, largely because of the greater average line 
width {crFWHM,W5NW ~ 3.5 km s~"^). As with the whole 
of W5, there is far too much turbulent energy in W5NW 
to be provided by outflows alone, implying the presence of 
another driver of turbulence. 

Figure [26] shows the W5NWpc region and Perseus 
molecular cloud on the same scale, though in two differ- 
ent emission lines. The Perseus cloud contains many more 
outflows and candidates (70 in Perseus vs. 13 in W5NWpc) 
despite a much larger physical area surveyed in W5. While 
it is likely that many of the W5W outflows will break apart 
into multiple flows at higher resolution, it does not seem 
likely that each would break apart into 5 flows, as would 
be required to bring the numbers into agreement. Since the 
highest density of outflows in Perseus is in the NGC 1333 
cloud, it may be that there is no equivalently evolved region 
in W5NWpc. The W5W region may be comparably mas- 
sive, but it is also confused and strongly interacting with 
the W5 HII region - either star formation is suppressed in 
this region, or outflows are rendered undetectable. In the 
latter case, the most likely mechanisms for hiding outflows 
are molecular dissociation by ionizing radiation and velocity 
confusion. 

Another possibility highlighted in Figure [26] is that the 
W5NW region is interacting with the W4 bubble. The cloud 
in the top right of Figure [26] is somewhat cometary, has 
higher peak brightness temperature, and is at a slightly dif- 



ferent velocity (-45 km s~^) than W5NW. The velocity dif- 
ference of 8 km s~^ could simply be two clouds physically 
unassociated along the line of sight, or could indicate the 
presence of another expanding bubble pushing two sheets 
of gas away from each other. Either way, the northwest por- 
tion of the W5NW region is clumpier than the area in which 
the outflows were detected, and it includes no outflow de- 
tections. 



5.2 Star Formation Activity 

CO outflows are an excellent tracer of ongoing embedded 
star formation (e.g.ll Shu et al.|1987 ). We use the locations of 
newly discovered outflows to qualitatively describe the star 
formation activity within the W5 complex and evaluate the 
hypothesis that star formation has been triggered on small 
or intermediate scales. 

Class O/I objects are nearly always associated with out- 
flows in nearby star- forming regions (e.g. Perseus Curtis 
Hatchell et al. 2007). However, Koenig et al. 



et al. 



2010 



(2008 ) detected 171 Class I sources in W5 using Spitzer pho- 



tometry. Since our detection threshold for outflow appears 



to be similar to that in Perseus (Section 3.1.1), the lower 
number of outflow detections is surprising, especially con- 
sidering that some of the detected outflows are outside the 
Spitzer-MIPS field (MIPS detections are required for Class 
I objects, and fiows 1-4 are outside that range) or are in 
the outer arm (flows 39-54). Additionally, we should detect 
outflows from Class objects that would not be identified 
by Spitzer colors. 

There are a number of explanations for our detection 
deficiency. The Class I objects detected within the HII re- 
gion "bubble" most likely have outfiows in which the CO is 
dissociated similar to jet systems in Orion (e.g. HH46/47, a 
pc-scale flow in which CO is only visible very near the pro- 
tostar; [Chernin k MassonlflOQll [Stanke et "al][1999t . This 
hypothesis can be tested by searching for optical and in- 
frared jets associated with these objects, which presumably 
have lower mass envelopes and therefore less extinction than 
typical Class I objects. Additionally, there are many outflow 
systems that are are likely to be associated with clusters of 
outflows rather than individual outflows as demonstrated in 
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(b) 

Figure 26. (a) An integrated CO 3-2 image of the W5W/NW region with ehipses overlaid displaying the locations and sizes of outflows. 
The dark red and blue ellipses in the lower right are associated with outer-arm outflows. W5W is the bottom- left, CO-bright region. 
W5NW is the top-center region containing the cluster of outflows, (b) An integrated CO 1-0 image of the Perseus molecular cloud from 
the COMPLETE survey ( |Arce et al.||201Q| . Note that the spatial scale is identical to that of (a) assuming that W5 is 8 times more 
distant than Perseus. The green ellipses represent |Arce et al.| ( 12010^ ) CPOCs while the orange represent known outflows from the same 
paper. 



Section [3.1.11 where we were able to identify fewer outflows 
when 'observing' the Perseus objects at a greater distance. 



There are 24 sources in the Koenig et al. (2008) Class I cat- 
alog within 15'' (one JCMT beam at 345 GHz) of another , 
and in many cases there are multiple Koenig et al. (2008) 



Class I sources within the contours of a single outflow sys- 
tem. 



5.3 Evaluating Triggering 

In the previous section, we discussed in detail the relation- 
ship between each sub-region and the HII region. Some re- 
gions are observed to be star-forming but not interacting 
with the HII region (W5NW, Sh 2-201), while others are 
interacting with the HII region but show no evidence or re- 
duced evidence of star formation (W5SE, W5W, LW Cas). 
At the very least, there is significant complexity in the trig- 
gering mechanisms, and no one mechanism or size scale is 
dominant. If we were to trust outflows as unbiased tracers 
of star formation, we might conclude that the majority of 



star formation in W5 is untriggered (spontaneous) , but such 
a conclusion is unreliable because both radiatively triggered 
star formation and "revealed" star formation may not ex- 
hibit molecular outflows (although ionized atomic outflows 
should still be visible around young stars formed through 
these scenarios). 

In Section |4.3| we analyzed a particular case in which 
the RDI mechanism could plausibly have crushed a cloud to 
create the observed protostar. It is not possible to determine 
whether interaction with the HII region was a necessary pre- 
condition for the star's formation, but it at least accelerated 
the process. The other cometary clouds share this property, 
but in total there are only 5 cometary clouds with detected 
outflows at their tips, indicating that this mechanism is not 
the dominant driver of star formation in W5. 

The 'collect and collapse' scenario might naively be ex- 
pected to produce an excess of young stars at the interac- 
tion front between the HII region and the molecular cloud. 
However, because such interactions naturally tend to form 
instabilities, this scenario produces cloud morphologies in- 
distinguishable from those of RDI. There is not an obvious 
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excess of sources associated with cloud edges over those deep 
within the clouds (e.g., Figure |2|. We therefore cannot pro- 
vide any direct evidence for this triggering scenario. 

The overall picture of W5 is of two concurrent episodes 
of massive- star formation that have lead to adjacent blown- 
out bubbles. Despite the added external pressure along the 
central ridge, it is relatively deficient in both star formation 
activity and dense gas, perhaps because of heating by the 
strong ionizing radiation field. The lack of star formation 
along that central ridge implies that much of the gas was 
squeezed and heated, but it was not crushed into gravita- 
tionally unstable fragments. While some star formation may 
have been triggered in W5, there is strong evidence for pre- 
existing star formation being at least a comparable, if not 
the dominant, mechanism of star formation in the complex. 



6 OUTFLOW SYSTEMS BEYOND W5 

Fifteen outflows were detected at velocities inconsistent with 
the local W5 cloud velocities. Of these, 8 are consistent with 
Perseus arm velocities (vlsr > — 55 km s~^) and could be 
associated with different clouds within the same spiral arm. 
The other 7 have central velocities vlsr < — 55 km s~^ and 
are associated with the outer arm identified in previous sur- 
Digel et al.|1996) . The properties of these outflows 



veys (e.g 

are given in Tables |4] and [5] the distances listed are kinematic 
distances assuming Rq = 8.4 kpc and vq = 254 km s~^ ([ReidJ 



et al.|[20Q9J . 

Of these outflows, all but one are within 2' of an IRAS 
point source. Outflow 54 is the most distant in our survey 
at a kinematic distance d = 7.5 kpc {visr — —75.6 km s""*^) 
and galactocentric distance Dq — 14.7 kpc. It has no known 
associations in the literature. 

Outflows 41 - 44 are associated with a cloud at vlsr ^ 
—62 km s~^ known in the literature as LDN 1375 and as- 
sociated with IRAS 02413+6037. Outflows 53 and 55 are 
at a similar velocity and associated with IRAS 02598+6008 
and IRAS 02425+6000 respectively. All of these sources he 
roughly on the periphery of the W5 complex. 

Outflows 45 - 52 are associated with a string of IRAS 
sources and HII regions to the north of W5 and have ve- 
locities in the range —55 < vlsr < —45. They therefore 
could be in the Perseus arm but are clearly unassociated 
with the W5 complex. Outflows 45 and 46 are associated 
with IRAS 02435+6144 and they may also be associated 
with the Sharpless HII region Sh 2-194. Outflows 47 and 48 
are associated with IRAS 02461+6147, also known as AFGL 
5085. Outflows 49 and 50 are nearby but not necessarily 
associated with IRAS 02475+6156, and may be associated 
with Sh 2-196. Outflows 51 and 52 are associated with IRAS 
02541+6208. 



7 CONCLUSIONS 

We have identified 40 molecular outflow candidates in the 
W5 star forming region and an additional 15 outflows spa- 
tially coincident but located in the outer arm of the Galaxy. 

• The majority of the CO clouds in the W5 complex are 
forming stars. Star formation is not limited to cloud edges 
around the HII region. Because star formation activity is 



observed outside of the region of influence of the W5 O- 
stars, it is apparent that direct triggering by massive star 
feedback is not responsible for all of the star formation in 
W5. 

• The W5 complex is seen nearly face-on as evidenced by 
a strict upper limit on the CO column through the center of 
the Hll-region bubbles. It is therefore an excellent region to 
study massive star feedback and revealed and triggered star 
formation. 

• Outflows contribute negligibly to the turbulent energy 
of molecular clouds in the W5 complex. This result is un- 
surprising near an HII region, but supports the idea that 
massive star forming regions are qualitatively different from 
low-mass star-forming regions in which the observed turbu- 
lence could be driven by outflow feedback. Even in regions 
far separated from the 0-stars, there is more turbulence and 
less energy injection from outflows than in, e.g., Perseus. 

• Despite detecting a significant number of powerful out- 
flows, the total outflowing mass detected in this survey 

1.5 Mq without optical depth correction, perhaps 10 — 20 
Mq when optical depth is considered) was somewhat smaller 
than in Perseus, a low to intermediate mass star forming re- 
gion with ~ 1/6 the molecular mass of W5. 

• The low mass measured may be partly because the CO 
3-2 line is sub-thermally excited in outflows. Therefore, while 
CO 3-2 is an excellent tracer of outflows for detection, it does 
not serve as a 'calorimeter' in the same capacity as CO 1-0. 

• Even considering excitation and optical depth correc- 
tions, it is likely that the mass of outflows in W5 is less than 
would be expected from a simple extrapolation from Perseus 
based on cloud mass. CO is likely to be photodissociated in 
the outflows when they reach the HII region, accounting 
for the deficiency around the HII region edges. However, in 
areas unaffected by the W5 0-stars such as W5NW, the de- 
ficiency may be because the greater turbulence in the W5 
clouds suppresses star formation or hides outflows. 

• Velocity gradients across the tails of many cometary 
clouds have been observed, hinting at their geometry and 
confirming that the outflows seen from their heads must be 
generated by protostars within. 

• Outflows have been detected in the Outer Arm at galac- 
tocentric distances > 12 kpc. These represent some of the 
highest galactocentric distance star forming regions discov- 
ered to date. 
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Table 2: CO 3-2 Outflow Measured Properties 



Outflow Latitude Longitude Ellipse Ellipse Ellipse Velocity Velocity J TXdv Bipolar? 

Number Major Minor PA center min max 

(kms-^) (kms-^) (kms-^) (K km s"^) 



lb 


136.4437 


1.2622 


60 


27 


342 


-36.1 


-47.6 


-40.3 


1.0 


yc 


Ir 


136.4674 


1.2705 


49 


24 


346 


-36.1 


-31.9 


-23.4 


1.5 


yc 


2b 


136.4899 


1.1904 


30 


23 


299 


-35.7 


-48.0 


-39.7 


0.7 


yc 


2r 


136.4743 


1.2042 


31 


28 


332 


-35.7 


-31.7 


-23.0 


1.3 


yc 


3 


136.475 


1.2548 


35 


25 


332 


-31.8 


-31.8 


-26.8 


1.3 


n 


4b 


136.5038 


1.2623 


26 


22 


35 


-36.2 


-44.1 


-40.1 


0.8 


yu 


4r 


136.5109 


1.2751 


25 


22 


332 


-36.2 


-32.4 


-28.6 


0.9 


yu 


5r 


136.5126 


1.2453 


24 


22 


10 


-35.3 


-31.4 


-28.8 


0.8 


yu 


5b 


136.5236 


1.2524 


39 


22 


3 


-35.3 


-45.0 


-39.2 


1.4 


yu 


6b 


136.532 


1.228 


28 


25 


332 


-35.3 


-44.8 


-40.0 


0.4 


yc 


6r 


136.5327 


1.2333 


28 


20 


318 


-35.3 


-30.6 


-24.0 


1.0 


yc 


7b 


136.5453 


1.2318 


24 


19 


332 


-34.9 


-47.5 


-39.9 


1.7 


yc 


7r 


136.5506 


1.2383 


27 


23 


314 


-34.9 


-29.9 


-22.7 


1.3 


yc 


8b 


136.5799 


1.2755 


18 


14 


332 


-34.5 


-41.5 


-39.3 


0.6 


yc 


8r 


136.581 


1.2601 


34 


30 


332 


-34.5 


-29.6 


-23.9 


1.4 


yc 


9b 


136.67 


1.2123 


30 


27 


332 


-35.0 


-44.5 


-38.5 


1.4 


yc 


9r 


136.6766 


1.2059 


40 


31 


332 


-35.0 


-31.6 


-26.7 


0.3 


yc 


10b 


136.7172 


0.7859 


39 


24 


353 


-42.8 


-52.6 


-47.5 


3.3 


yc 


lOr 


136.7271 


0.7797 


31 


26 


332 


-42.8 


-38.1 


-33.1 


4.1 


yc 


lib 


136.8195 


1.082 


25 


24 


331 


-34.2 


-40.7 


-37.0 


3.1 


yc 


llr 


136.8173 


1.0799 


24 


22 


331 


-34.2 


-31.4 


-20.4 


1.5 


yc 


12b 


136.8414 


1.1512 


30 


26 


332 


-40.4 


-53.3 


-46.2 


1.5 


yc 


12r 


136.8479 


1.1517 


27 


25 


332 


-40.4 


-34.6 


-30.1 


0.9 


yc 


13 


136.8461 


0.8426 


28 


27 


332 


-31.0 


-31.0 


-23.5 


1.0 


n 


14 


136.8591 


1.176 


24 


23 


332 


-47.1 


-54.5 


-47.1 


0.8 


n 


15 


136.9443 


1.0841 


28 


18 


348 


-45.0 


-55.0 


-45.0 


3.1 


n 


16b 


137.3929 


0.5977 


23 


18 


333 


-40.7 


-47.0 


-42.6 


0.7 


yu 


16r 


137.3981 


0.6121 


22 


19 


357 


-40.7 


-38.7 


-35.2 


1.9 


yu 


17b 


137.4084 


0.6762 


20 


18 


293 


-40.3 


-57.9 


-43.0 


2.3 


yc 


17r 


137.412 


0.6775 


20 


18 


308 


-40.3 


-37.6 


-30.4 


1.1 


yc 


18b 


137.4925 


0.6289 


16 


15 


333 


-35.5 


-39.2 


-37.6 


1.1 


yc 


18r 


137.4908 


0.6292 


18 


17 


307 


-35.5 


-33.4 


-31.0 


2.0 


yc 


19b 


137.4815 


0.6409 


20 


17 


1 


-36.0 


-41.9 


-38.9 


1.3 


yc 


19r 


137.4798 


0.6404 


20 


16 


301 


-36.0 


-33.1 


-25.9 


0.7 


yc 


20r 


137.5368 


1.2792 


24 


21 


332 


-37.4 


-33.0 


-22.5 


5.2 


yc 


20b 


137.539 


1.279 


27 


23 


17 


-37.4 


-52.0 


-41.8 


3.4 


yc 


21b 


137.6152 


1.3543 


31 


28 


322 


-39.5 


-52.0 


-43.7 


4.5 


yc 


21r 


137.6169 


1.3585 


31 


18 


4 


-39.5 


-35.2 


-30.0 


1.2 


yc 


22 


137.6213 


1.506 


27 


21 


293 


-40.3 


-46.0 


-40.3 


2.1 


n 


23b 


137.6389 


1.5251 


21 


14 


331 


-38.5 


-42.5 


-40.5 


1.6 


yc 


23r 


137.6449 


1.5194 


19 


12 


331 


-38.5 


-36.5 


-32.0 


1.9 


yc 


24r 


137.7094 


1.4824 


20 


20 


331 


-38.2 


-33.8 


-25.4 


4.2 


yc 


24b 


137.7146 


1.4809 


25 


19 


292 


-38.2 


-50.0 


-42.7 


4.4 


yc 


25b 


138.1398 


1.6858 


39 


26 


282 


-38.8 


-49.5 


-43.2 


0.6 


yc 


25r 


138.142 


1.6884 


43 


35 


11 


-38.8 


-34.3 


-27.5 


1.7 


yc 


26b 


138.2913 


1.5538 


29 


29 


355 


-38.7 


-52.0 


-47.4 


1.2 


yc 


26r 


138.2966 


1.5564 


28 


28 


330 


-38.7 


-30.0 


-20.0 


4.2 


yc 


27 


138.3017 


1.5689 


26 


25 


330 


-30.0 


-30.0 


-22.0 


1.8 


n 


28 


138.3042 


1.5437 


20 


19 


330 


-43.3 


-46.1 


-43.3 


1.4 


n 


29 


138.3053 


1.5537 


22 


20 


330 


-45.3 


-51.6 


-45.3 


2.5 


n 


30 


138.3115 


1.5443 


26 


26 


330 


-33.0 


-33.0 


-29.2 


1.2 


n 


31 


138.3184 


1.5566 


26 


25 


330 


-44.4 


-49.1 


-44.4 


1.1 


n 



Continued on next page 
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Outflow 


Latitude 


Longitude 


Ellipse 


Ellipse 


Ellipse 


Velocity 


Velocity 




fTXdv 


Bipolar? ^ 


Number 






Major 


Minor 


PA 


center 


min 


max 












II 


// 


o 


(km s-^) 


(km s-^) 


(km s-^) 


(K km s-^) 




32 


138.3213 


1.5658 


27 


27 


330 


-31.7 


-31.7 


-27.0 


1.4 


n 


33b 


138.3618 


1.5073 


28 


26 


330 


-39.4 


-49.5 


-44.0 


1.3 


yc 


33r 


138.3642 


1.4959 


29 


21 


330 


-39.4 


-34.7 


-25.8 


2.0 


yc 


34r 


138.4779 


1.6137 


22 


21 


330 


-36.9 


-33.1 


-29.1 


0.5 


yc 


34b 


138.4768 


1.6142 


21 


20 


330 


-36.9 


-43.6 


-40.6 


0.8 


yc 


35r 


138.4998 


1.6496 


22 


20 


4 


-37.5 


-31.3 


-24.1 


1.4 


yc 


35b 


138.5021 


1.6458 


23 


21 


330 


-37.5 


-49.5 


-43.6 


1.3 


yc 


36b 


138.5034 


1.6654 


35 


26 


5 


-37.5 


-50.4 


-42.3 


1.2 


yc 


36r 


138.5061 


1.6576 


22 


21 


330 


-37.5 


-32.6 


-26.7 


1.4 


yc 


37r 


138.5208 


1.6618 


27 


22 


330 


-38.5 


-33.5 


-31.4 


0.6 


yc 


37b 


138.5241 


1.6667 


23 


23 


18 


-38.5 


-47.0 


-43.6 


0.6 


yc 


38b 


137.4983 


0.6062 


16 


15 


333 


-36.1 


-39.2 


-38.5 


0.8 


yc 


38r 


137.4977 


0.6055 


15 


14 


307 


-36.1 


-33.7 


-32.5 


0.5 


yc 


39b 


138.1506 


0.7724 


23 


16 


321 


-38.8 


-45.3 


-41.0 


2.0 


yc 


39r 


138.1591 


0.7713 


17 


13 


304 


-38.8 


-36.6 


-34.7 


0.7 


yc 


40 


138.1356 


0.7634 


22 


18 


4 


-36.0 


-36.0 


-27.6 


2.2 


n 



Measured properties of the outflows. 

" Is the outflow part of a bipolar pair? yc = yes, confident; yu = yes, uncertain; n = no 
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Table 3: CO 3-2 Outflow Derived Properties 



Outflow Mass Momentum Energy Dynamical Momentum 
Number Age Flux 

(Mo) (Mokm s-^) (10^2 ergs) (10^ years) 10"^ M© 

km s~^yr~^ 



lb 


0.034 


0.26 


21.1 


7.0 


7.2 


Ir 


0.04 


0.24 


17.3 


7.0 


7.2 


2b 


0.011 


0.07 


4.9 


5.4 


4.4 


2r 


0.025 


0.17 


13.0 


5.4 


4.4 


3 


0.025 


0.12 


5.8 


_ 


_ 


4b 


0.01 


0.06 


4.0 


7.2 


1.5 


4r 


0.011 


0.04 


1.8 


7.2 


1.5 


5r 


0.01 


0.04 


2.0 


4.5 


4.0 


5b 


0.025 


0.14 


8.0 


4.5 


4.0 


6b 


0.007 


0.04 


3.0 


1.7 


8.1 


6r 


0.013 


0.09 


6.8 


1.7 


8.1 


7b 


0.017 


0.13 


10.5 


2.4 


10.9 


7r 


0.018 


0.13 


9.7 


2.4 


10.9 


8b 


0.003 


0.02 


0.9 


4.9 


4.5 


8r 


0.032 


0.2 


13.2 


4.9 


4.5 


9b 


0.025 


0.13 


7.2 


3.9 


4.2 


9r 


0.009 


0.04 


1.8 


3.9 


4.2 


10b 


0.068 


0.41 


25.7 


3.9 


22.0 


lOr 


0.074 


0.45 


28.0 


3.9 


22.0 


lib 


0.042 


0.17 


7.2 


0.7 


35.3 


llr 


0.017 


0.09 


5.9 


0.7 


35.3 


12b 


0.026 


0.14 


8.7 


1.8 


15.2 


12r 


0.014 


0.13 


12.7 


1.8 


15.2 


13 


0.016 


0.1 


5.8 


_ 


_ 


14 


0.01 


0.06 


4.1 


_ 


_ 


15 


0.036 


0.24 


17.3 


_ 


_ 


16b 


0.006 


0.03 


1.2 


11.1 


0.7 


16r 


0.018 


0.05 


1.3 


11.1 


0.7 


17b 


0.019 


0.12 


9.4 


1.4 


10.6 


17r 


0.009 


0.03 


0.7 


1.4 


10.6 


18b 


0.006 


0.02 


0.5 


1.4 


4.0 


18r 


0.013 


0.04 


1.0 


1.4 


4.0 


19b 


0.011 


0.05 


2.4 


0.7 


9.6 


19r 


0.005 


0.01 


0.2 


0.7 


9.6 


20r 


0.059 


0.5 


46.3 


0.5 


156.0 


20b 


0.047 


0.33 


26.6 


0.5 


156.0 


21b 


0.086 


0.58 


41.4 


1.7 


39.1 


21r 


0.014 


0.08 


4.3 


1.7 


39.1 


22 


0.027 


0.1 


4.3 






23b 


0.011 


0.03 


0.9 


4.5 


1.3 


23r 


0.01 


0.03 


1.0 


4.5 


1.3 


24r 


0.037 


0.3 


26.1 


1.7 


34.1 


24b 


0.047 


0.28 


18.3 


1.7 


34.1 


25b 


0.014 


0.09 


6.8 


1.0 


42.8 


25r 


0.056 


0.35 


23.0 


1.0 


42.8 


26b 


0.023 


0.24 


26.1 


1.1 


98.3 


26r 


0.072 


0.85 


106.0 


1.1 


98.3 


27 


0.026 


0.09 


4.5 






28 


0.012 


0.07 


4.6 






29 


0.024 


0.06 


2.1 






30 


0.018 


0.12 


8.0 







Continued on next page 
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Outflow Mass Momentum 
Number 

(Mo) (Mokm s-^) 



Energy Dynamical 



(10^^ ergs) (10^ years) 



Momentum 
Flux 

10-^ Mo 
km s~"'^yr~"'^ 



31 


0.016 


0.03 


0.7 






32 


0.023 


0.18 


14.5 






33b 


0.022 


0.14 


10.1 


3.0 


11.4 


33r 


0.026 


0.2 


16.1 


3.0 


11.4 


34r 


0.005 


0.03 


1.7 


0.7 


8.4 


34b 


0.007 


0.03 


1.2 


0.7 


8.4 


35r 


0.013 


0.12 


11.6 


1.3 


18.7 


35b 


0.014 


0.12 


11.0 


1.3 


18.7 


36b 


0.025 


0.19 


15.8 


2.7 


10.7 


36r 


0.014 


0.1 


6.8 


2.7 


10.7 


37r 


0.008 


0.04 


1.6 


2.1 


4.3 


37b 


0.007 


0.06 


4.2 


2.1 


4.3 


38b 


0.005 


0.01 


0.4 


0.9 


2.3 


38r 


0.002 


0.01 


0.1 


0.9 


2.3 


39b 


0.017 


0.07 


2.8 


7.5 


1.0 


39r 


0.004 


0.01 


0.3 


7.5 


1.0 


40 


0.019 


0.08 


3.5 







Derived properties of the outflows in the optically thin limit. 

Typical optical depth corrections for ^^C03-2 are in the range 7-14 jCurtis et al. 20101. 

The correction for velocity confusion is > 2 but poorly constrained (Arce et al. 12010). 

Finally, an excitation correction in the range 1-20 is likely required as described in tne Appendix. 

The mass and momentum values can be multiplied by these factors to acquire the corrected values. 

The energy is weighted more heavily towards high- velocity, low-optical-depth gas, so the correction factor is likely to be lower. 
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Table 4. Outer Arm CO 3-2 Outflows - Measured Properties 



1 1 f fl CWM 
Ky Ll Liii^ vv 




T ,0"n ci'l'iirip 


Rllinse 


Rllinse 


Rllinse 




IXQ 


V C-i^ v^l L V 






J ->- 


Number 






]V[£ijor 


IVIinor 


PA 


Distance 


















II 


II 


o 






(km s~^) 






(K km s^l") 


41r 


136.364 


0.9606 


25 


18 


2 


5510 


13000 


-61.8 


-59.2 


-56.5 


0.5 


41b 


136.3634 


0.9568 


23 


17 


353 


5510 


13000 


-61.8 


-71.6 


-64.3 


3.0 


42r 


136.3522 


0.9786 


20 


14 


2 


5500 


12900 


-62.1 


-59.8 


-57.6 


0.6 


42b 


136.3548 


0.9798 


20 


19 


332 


5500 


12900 


-62.1 


-67.8 


-64.4 


0.5 


43r 


136.3495 


0.9612 


17 


15 


63 


5510 


13000 


-61.8 


-59.0 


-56.1 


0.8 


43b 


136.353 


0.9621 


12 


12 


333 


5510 


13000 


-61.8 


-66.3 


-64.6 


1.0 


44r 


136.3554 


0.9576 


13 


13 


23 


5500 


12900 


-61.8 


-59.0 


-55.4 


2.1 


44b 


136.3545 


0.9567 


14 


14 


333 


5500 


12900 


-61.8 


-68.0 


-64.5 


2.0 


45r 


136.1219 


2.0816 


34 


25 


297 


3750 


11400 


-46.5 


-43.1 


-40.5 


0.6 


45b 


136.1233 


2.0803 


35 


25 


306 


3750 


11400 


-46.5 


-57.3 


-50.0 


1.9 


46 


136.1166 


2.0983 


26 


25 


332 


3790 


11400 


-50.2 


-52.6 


-50.2 


0.5 


47b 


136.3857 


2.2687 


34 


27 


332 


3220 


11000 


-42.0 


-55.0 


-46.7 


3.5 


47r 


136.3861 


2.267 


35 


23 


304 


3220 


11000 


-42.0 


-37.3 


-25.1 


5.0 


48b 


136.374 


2.2628 


29 


21 


332 


3250 


11000 


-43.2 


-51.4 


-47.0 


1.5 


48r 


136.3736 


2.2615 


29 


22 


332 


3250 


11000 


-43.2 


-39.5 


-22.2 


8.9 


49r 


136.4663 


2.4678 


29 


23 


290 


3610 


11300 


-45.7 


-42.2 


-33.0 


2.2 


49b 


136.4661 


2.4693 


31 


23 


292 


3610 


11300 


-45.7 


-52.2 


-49.1 


0.9 


50b 


136.5087 


2.5108 


31 


25 


332 


3380 


11100 


-43.5 


-48.5 


-46.5 


0.8 


50r 


136.5118 


2.5083 


28 


23 


10 


3380 


11100 


-43.5 


-40.6 


-37.5 


1.0 


51b 


137.058 


2.9858 


28 


23 


293 


4350 


11900 


-51.8 


-55.5 


-53.0 


0.8 


51r 


137.0567 


2.9864 


34 


25 


8 


4350 


11900 


-51.8 


-50.6 


-40.9 


3.5 


52r 


137.0662 


2.9999 


37 


26 


43 


4390 


12000 


-52.2 


-49.1 


-41.0 


7.8 


52b 


137.0683 


3.0013 


38 


29 


15 


4390 


12000 


-52.2 


-65.8 


-55.2 


4.1 


53b 


138.6143 


1.5611 


26 


26 


330 


5450 


13000 


-59.7 


-71.1 


-61.7 


5.5 


53r 


138.6158 


1.563 


25 


23 


330 


5450 


13000 


-59.7 


-57.6 


-54.5 


1.3 


54r 


136.382 


0.8392 


29 


20 


343 


7480 


14700 


-75.6 


-73.2 


-68.9 


2.6 


54b 


136.3824 


0.838 


20 


17 


332 


7480 


14700 


-75.6 


-83.1 


-77.9 


2.0 


55b 


136.7623 


0.4548 


27 


16 


343 


5230 


12700 


-60.9 


-65.2 


-62.7 


5.4 


55r 


136.7579 


0.4522 


24 


18 


343 


5230 


12700 


-60.9 


-59.0 


-53.2 


6.0 



El 

"Galactocentric Radius 
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Table 5. Outer Arm CO 3-2 Outflows - Derived Properties 



Outflow 


Mass 


Alomentum 




Dyn£irnic£il 


JVLomentum 


Nuniber 








Affe 


Flux 




V^-'-Q ) 


(M0km s"^) 




(10^ yecirs) 


10-6 M0km s"^yr"^ 


41r 


0.037 


0.11 


3.5 


3.6 


30.2 


41b 


0.196 


0.96 


54.6 


3.6 


30.2 


42r 


0.029 


0.06 


1.3 


4.3 


3.9 


42b 


0.03 


0.11 


3.9 


4.3 


3.9 


43r 


0.033 


0.13 


5.1 


7.3 


2.9 


43b 


0.024 


0.08 


2.8 


7.3 


2.9 


44r 


0.062 


0.21 


7.9 


1.8 


27.7 


44b 


0.067 


0.28 


12.6 


1.8 


27.7 


45r 


0.037 


0.18 


8.6 


1.2 


62.3 


45b 


0.126 


0.55 


28.3 


1.2 


62.3 


46 


0.028 


0.1 


3.7 


- 


- 


47b 


0.187 


1.32 


101.0 


0.6 


553.0 


47r 


0.232 


1.84 


164.0 


0.6 


553.0 


48b 


0.054 


0.33 


20.8 


0.4 


754.0 


48r 


0.341 


2.4 


229.0 


0.4 


754.0 


49r 


0.106 


0.77 


62.7 


1.4 


71.1 


49b 


0.047 


0.22 


10.8 


1.4 


71.1 


50b 


0.037 


0.14 


5.5 


3.8 


7.4 


50r 


0.038 


0.14 


5.5 


3.8 


7.4 


51b 


0.058 


0.13 


3.0 


1.0 


152.0 


51r 


0.303 


1.33 


72.7 


1.0 


152.0 


52r 


0.829 


4.3 


250.0 


1.4 


479.0 


52b 


0.472 


2.5 


150.0 


1.4 


479.0 


53b 


0.626 


3.16 


194.0 


4.7 


73.2 


53r 


0.124 


0.31 


8.5 


4.7 


73.2 


54r 


0.461 


1.24 


37.5 


1.8 


135.0 


54b 


0.212 


1.14 


64.0 


1.8 


135.0 


55b 


0.411 


1.66 


68.0 


7.5 


28.3 


55r 


0.404 


0.47 


6.6 


7.5 


28.3 
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APPENDIX A: OPTICALLY THIN, LTE DIPOLE MOLECULE 

While many authors have solved the problem of converting CO 1-0 beam temperatures to H2 column densities (jGarden et al. 



1991||Bourke et al.|1997||Cabrit BertoutllQQOllLada k Fich|1996 |), there are no examples in the literature of a full derivation 



of the LTE, optically thin CO-to-H2 conversion process for higher rotational states. We present the full derivation here, and 
quantify the systematic errors generated by various assumptions. 

We begin with the assumption of an optically thin cloud such that the radiative transfer equation ( Wilson et al.|[2009 
eqn 1.9) simplifies to 



ds 

The absorption and stimulated emission terms yield 

h^yulBuln^^ 



= -Kulu (Al) 



hVulBluTll 

zy A2 

c c 

where ^{f) is the line shape function (J Lp{v)di' = 1), n is the density in the given state, v is the frequency of the transition, 
B is the Einstein B coefficient, and h is Planck's constant. 

By assuming LTE (the Boltzmann distribution) and using Kirchoff's Law and the definition of the Einstein A and B 
values, we can derive a more useful version of this equation 



exp 



1 



where /cs is Boltzmann's constant. 

The observable Tb can be related to the optical depth, which is given by 



exp 



Uuds 



(A3) 



(A4) 



Rearranging and converting from density to column (J nds = N) gives an equation for the column density of the molecule 
in the upper energy state of the transition: 



Nu 



c^Aui 



exp 



(A5) 



In order to relate the brightness temperature to the optical depth, at CO transition frequencies the full blackbody formula 
must be used and the CMB must also be taken into account. Wilson et al. ( 2009| equation 15.29 



TbW) = 



is rearranged to solve for Tu'. 



-In 



kB 



1 



ksTB ( \ hv 



"1 r hv 



/kuTcMB 



.]-')(. 



l]-'-[e- 



/ksTcMB 



We convert from frequency to velocity units with du = ly/cdv, and plug (A7) into (A5) to get 

-1 



exp 



ksTex 



In 



ksTB 



1 



/kBTcMB 



dv 



which is the full LTE upper- level column density with no approximations applied. 
The first term of the Taylor expansion is appropriate for r << 1 (ln[l + x] ^ x - 



Nu 



c^Aui 



exp 



\kBTe: 



J huui V L 



l^ul/ksTe: 



dv 



which simplifies to 



Nu 



STTulikE 



C^Aulh e^^ul/kBTcMB — ehiy^i/kBTe. 



I 



Tadv 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



This can be converted to use /ie (0.1222 for ^^CO; Muenter||l975| , the electric dipole moment of the molecule, instead 



of Aui , using 



Wilson et al. 



(2009) equation 15.20 {{Aui = (647r^)/(3/ic^)) ly^fil): 



Nu 



3 kB 2Ju + 1 



phiyul/kBTcmb 



1 



/ 



Tgdv 



(All) 



Stt^ Uui Ju c^'^^i/^bTcmb - e^^'^^l/^BTex 
The total column can be derived from the column in the upper state using the partition function and the Boltzmann 



distribution 



ntot 



nj — no 

J=0 J=0 



^(2J+l)exp( — — 



(A12) 
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10 20 30 40 50 60 

Excitation Temperature (K) 

Figure Al. The LTE, optically thin conversion factor from Tb (K km s-^) to N(H2) (cm-^) 

assuming Xi2qq = 10 ^ plotted against 

Tex- The dashed line shows the effect of using the integral approximation of the partition function (e.g. |Cabrit Bertout|1990| >. It is a 
better approximation away from the critical point, and is a better approximation for higher transitions. The dotted line shows the effects 
of removing the CMB term from ( |A6| ; the CMB populates the lowest two excited states, but contributes nearly nothing to the J = 3 
state. Top (blue): J=l-0, Middle (green): J=2-l, Bottom (red): J=3-2. 



This equation is frequently approximated using an integral (e.g. jCabrit & Bertout 1990), but a more accurate numerical 
solution using up to thousands of rotational states is easily computed 

j{j + l)Beh' 



nj 



J2 (2j + l)exp 



(2J+ l)exp 



J{J+l)Beh 

ksTex 



(A13) 



The effects of using the approximation and the full numerical solution are shown in figure [Al] 

The CO 3-2 transition is also less likely to be in LTE than the 1-0 transition. The crit ical density (ricr = Aui/Cui) of 
^^CO 3-2 is 27 times higher than that for 1-0. We have run RADEX ( [van der Tak et al.|2007| LVG models of CO to examine 
the impact of sub-thermal excitation on column derivation. The results of the RADEX models are shown in Figure [A2] They 
illustrate that, while it is quite safe to assume the CO 1-0 transition is in LTE in most circumstances, a similar assumption 
is probably invalid for the CO 3-2 transition in typical molecular cloud environments. 
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Figure A2. Top: The derived N(H2) as a function of for Tb = 1 K. The dashed hnes represent the LTE-derived N{Yi2) /Tb factor, 
which has no density dependence and, for CO 3-2, only a weak dependence on temperature. We assume an abundance of ■'^^CO relative 
to H2 Xco — lO"'^. Bottom: The correction factor (N{irL2)RADEX / N(H2)lt£;) as a function of nu^. For Tk = 20 K, the "correction 
factor" at 10^ cm~^ (typical GMC mean volume densities) is ~ 15, while at 10^ cm""^ (closer to ricrit but perhaps substantially higher 
than GMC densities) it becomes negligible. The correction factor is also systematically lower for a higher gas kinetic temperature. For 
some densities, the "correction factor" dips below 1, particularly for CO 1-0. This effect is from a slight population inversion due to fast 
spontaneous decay rates from the higher levels and has been noted before (e.g. |Goldsmith|1972|>. 



